Asymptotic behavior of self-affine processes in semi- infinite domains 
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We propose to model the stochastic dynamics of a polymer passing through a pore (translocation) 
by means of a fractional Brownian motion, and study its behavior in presence of an absorbing 
boundary. Based on scaling arguments and numerical simulations, we present a conjecture that 
provides a link between the persistence exponent 6 and the Hurst exponent H of the process, thus 
sheding light on the spatial and temporal features of translocation. Furthermore, we show that 
this conjecture applies more generally to a broad class of self affine processes undergoing anomalous 
diffusion in bounded domains, and we discuss some significant examples. 



The dynamics of a polymer chain composed of N 
monomers passing through a pore (translocation) has 
been intensively investigated in recent years, by virtue 
of its central role in understanding, e.g., viral injection 
of DNA into a host or RNA transport through nanopores, 
and mastering such applications as fast DNA or RNA se- 
quencing through engineered channels 
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translocation coordinate s(t), namely the label of the 
monomer crossing the pore at time t, is key to quan- 
titatively describing the translocation process 0, [y, 0], 
which begins when s = 1 and ends when s = N, i.e., 
when the first and the last monomer of the chain enter 
the pore, respectively (see Fig. [T]). 

Various dynamical regimes of s(t) have been identi- 
fied: in absence of driving forces and hydrodynamic ef- 
fects (free polymer), fluctuations dominate and s(t) can 
be regarded as a stochastic process, whose features vary 
with the polymer length N [6|, [8|, l2( ■ Understanding the 
dynamics of s(t) represents a challenging problem. A 
free polymer is characterized by two natural time scales. 
First, the intrinsic equilibration time t eq required by the 
center of mass of the polymer R cm to travel a distance of 
the order of the typical size of the chain. This size is given 
by the radius of gyration R g , which scales as R g ~ N v 
in the large N limit. In a good solvent, v = 3/4 in 2d 
and v ~ 0.59 in 3d when exluded-volume effects for the 
monomers are considered; v = 1/2 for an ideal ('phan- 
tom') polymer [l(|. The center of mass diffuses with a 
diffusion coefficient ~ 1/N. Then, R cm ~ y/t/N . Hence, 
y/teg/N ~ Rg and t eq ~ N 2y+1 for large N. On the other 
hand, the translocation time T (much longer than t eq ) is 
the time required by the polymer to go through the pore, 
so that s(T) — N (Fig. [TJ. Under the hypothesis that 
the translocation is a self-affine process, i.e., s(t) ~ t H , 
with Hurst exponent H, it follows that T ~ N x / H . 

For short polymers, excluded- volume effects are neg- 
ligible, s(t) undergoes diffusion, and T ~ A^ 2 How- 
ever, as N increases, the excluded-volume interactions 
become relevant and s(t) undergoes subdiffusion, < 
H < 1/2 H. Numerical simulations (mostly 2d) support 
the following conclusions: i) T and t eq have the same 
scaling, up to a large prefactor, i.e., T ~ t eq . Hence, 




FIG. 1: Left: Translocation of a polymer chain through a 
pore. Right: The translocation coordinate s(t) denotes the 
number of the monomer that is crossing the pore at time t. 

H = 1/(1 + 2v) [8(; (Note that for 'phantom' polymers 
s(t) always diffuses, even for large N.) ii) the probabil- 
ity P(s, t) of finding the monomer s in the pore at time 
t for an infinite chain (i.e., in absence of boundaries) is 
Gaussian [111 ] ; Hi) for a finite chain yet to have completed 
translocation, s(t) evolves in presence of two absorbing 
boundaries at s = and s = N (see Fig. [T]) and the 
distribution of s(t) converges to a non-Gaussian form at 
long times. In particular, this distribution vanishes non- 
linearly as and (N — s)^ at the two boundaries s = 
and s = N respectively, with <f) ~ 1.44 in 2d [ll[. Com- 
puting analytically the exponent 0is a nontrivial problem 
and is the main objective of this Letter. 

Upon gathering these hints from simulations, we pro- 
pose as a natural candidate for s(t) the fractional Brow- 
nian motion (fBm), a self-affine Gaussian process (0 < 
H < 1). A Gaussian process is completely defined by its 
autocorrelation function (s(ti)s(t2)): for fBm, [l2| 

(s(h)s(t 2 )) <xtl H + t 2 2 H -\h-t 2 \ 2H . (1) 

Brackets refer to ensemble average over many realiza- 
tions. Equation [1] implies that the incremental correla- 
tion function of fBm is stationary 

((s(*i)-^2)) 2 )oc|^-i 2 | 2 ^. (2) 

Note that Brownian motion (BM) is fBm with H = 1/2. 
When H = 1/(1 + 2v), fBm satisfies conditions i) and 
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ii). Concerning Hi), the central result of our paper is to 
show that the probability distribution of a fBm confined 
in the positive half-axis vanishes as s^, with 
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close to the absorbing boundary s = 0, in the long time 
limit. For 2d polymers with excluded-volume effects, us- 
ing v — 3/4 one gets H = 2/5 and 4> = 3/2, which is in 



good agreement with the numerical value (j) ~ 1.44 11|. 
In 3d, using v ~ 0.59, we predict H ~ 0.46 and <f> ~ 1.18. 

More generally, we will show that actually for a broad 
class of self-affine processes with Hurst exponent H the 
exponent <f> is related to the persistence exponent 6 of the 
process via the relation tfi — 8/H. The persistence of a 
stochastic process is simply the probability of no return 
to its initial value up to time t and for a wide class of 
processes the persistence decays algebraically ~ t~ 6 at 



late times witha nontrivial persistence exponent 9 13 1 . 
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For fBm, in particular, 8 = 1 — H is known exactly 
leading to the result <\> = (1—H) / H . This general relation 
of <p to persistence exponent 9 sheds light on the spatial 
and temporal features of anomalous diffusion in presence 
of absorbing boundaries, which is ubiquitous in nature 
and arises is such diverse fields as contaminant migration 
in heterogeneous materials [15j and charge transport in 
amorphous semiconductors, to name only a few [161 ]. 

Let us begin by introducing the class of stochastic pro- 
cesses we are interested in. A generic unbiased stochastic 
walk x{t) can be represented as a sum of jumps (incre- 
ments) £(t'), 



(4) 



where the jump length has a symmetric marginal 
distribution 7iv(£). The Central Limit Theorem guaran- 
tees that, if x(t) undergoes anomalous diffusion, one of 
the following statements must be violated: a) the jumps 
are identically distributed, i.e., 7iy(£) does not depend on 
time: x(t) has stationary increments; b) the variables £ 
are independent: x(t) is Markovian; c) 7iv(£) has finite 
variance. For instance, fBm, which is the only self-affine 
Gaussian process with stationary increments, is anoma- 
lous and violates b) (except for H = 1/2). In the follow- 
ing, we present a general scaling argument for self-affine 
processes satisfying property a), i.e., having stationary 
increments, but not restricted to satisfy either b) or c). 

For such processes started in x$ > and killed upon 
leaving the positive half-axis, we define P+(x, x$, t) as 
the probability of finding the walker in x at time t. We 
define p Xo (x, t) as the conditional probability density of 
finding the walker, given that it has not been absorbed 
at any previous time: 



Px (x,t) 



P + (x,x ,t) 
J °° dxP + (x,x ,t) 



(5) 



At long times, the small- a; behavior of this distribution 
gives access to the exponent <f>. The quantity S(xo,t) — 
J °° dxP(x, xo,t) defines the survival probability that the 
walker has not left the positive half-axis up to t. 

As a useful guide, let us first recall the results for reg- 
ular BM, i.e., fBm with H = 1/2. The method of images 
gives the scaling form 17 1 
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where F(y,y ) = [e-fo-w) /2 _ e -(y+yo) /a 
long time behavior of the conditional probability, for any 
xq, is given by xq — > 0: 



P0 (x,i) = ^= e -* 2 / 2t . 



(7) 



Observe that in this limit the process is not reminiscent 
of the initial condition. For large x, the particles do not 
feel the presence of the boundary, and po(x, t) behaves as 
the Gaussian probability of an unconstrained BM. On the 
other hand, po(x,t) vanishes linearly close to the bound- 
ary, which implies that 0=1 for BM. 

Inspired by the special case of BM, we anticipate a 
similar scaling for a generic self-affine process: 



P+(x,x ,t) 
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where H' has to be determined. Integrating over x, 
we get S(x ,t) = t H - H 'G(x /t H ), where G(y) = 
J °° F(x, y)dx is a scaling function of a single variable. 
For any Xq > 0, we know that S(xq,0) = 1. Therefore, 
the only acceptable choice is H' = H , which implies 



S(x ,t) = G[^) 



(9) 



with G(+oo) = 1. This scaling form allows the asymp- 
totic behavior of S(xo,t) to be fully characterized. In 
particular, since the survival probability is expected to 
decay as S(xo,t) ~ t~ e for large t (the semi-infinite do- 
main lacks a characteristic scale), G(y) must behave as 
G(y) r-j y e / H for small y, so that 

S(x ,t)ott- e x e /H . (10) 

Consider now the scaled variables y = x/t and yo = 
xo/t H . By combining Eq. [9] with Eqs.[8]and[5l we get 

F(y,y )=S(y ) Pyo (y) (11) 

where p yo (y) is the conditional probability density ex- 
pressed in terms of the rescaled variables. In the long 
time limit, yo — > and F(y,yo) can be factorized as 



F {y,yo) oc vl /H po(y)- 



(12) 



Let us now consider the limit y — > and suppose that 
Po(y) ~ y? '. The process is time reversible since its incre- 
ments are stationary, i.e., a path from xq to x forward in 
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FIG. 2: Survival probability as a function of the rescaled vari- 
able j/o, for H — 4/7 (top) and H — 4/9 (bottom), at various 
times t in a semi-infinite domain. Numerical simulations are 
compared with the proposed scaling (1 — H)/H = 3/4 and 
(1 - H)/H = 5/4 respectively. 



time plays the same role of a path from x to xo backward 
in time. As a consequence, the two limits xq — > and 
x — ► can be interchanged in P+(x, xo, t), implying that 
yo and y must appear in an identical fashion in the lim- 
iting behavior of F(y,yo) in Eq. (TT2"|) . and thus follows 
our proposed scaling relation <p = 9/H. This conjecture 
represents our central result. 

In support of this conjecture, we proceed next to nu- 
merically confirm its validity for fBm. In order to gen- 
erate a fBm path of L + 1 steps, {x$, X\, Xi, xl}, 
we need to draw L Gaussian numbers correlated through 
Eq. [T] This is in general a time consuming procedure, as 
it requires computing the square root of the autocorre- 
lation matrix. Better results can be obtained by making 
use of the stationarity of the increments & = Xi — Xi-i, 
with a Toeplitz matrix correlator (£i£i+fc) cx \k — 1\ 2H + 
\k + l\ 2H — 2\k\ 2H . In particular, the Levinson algorithm 
allows first passage problems to be efficiently tackled, by 
recursively generating Xi+\ given Xi, . . . ,xo fill ]. 

In Fig. [2j we test the scaling of S(xo,t) proposed in 
Eq. 1101 The survival probability is shown as a function 
of the rescaled variable yo, for different times. Particles 
start at Xo > and a single absorbing boundary is set 
in x — 0. All curves collapse, in agreement with the self- 
affinity of the process. The collapse does not work in 
the region where xq is of the same order of magnitude as 
the typical increment of the (discrete) process. This ef- 
fect disappears as time (and thus path length) increases, 
approaching the continuum limit. A good agreement is 
found between the proposed conjecture and the numeri- 
cal results. In the superdiffusive regime, t H grows faster, 
and shorter times are required to observe the predicted 
scaling, i.e., (f> = (1 — H)/H. 
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FIG. 3: Semi- infinite domain. Conditional probability po(y) 
as a function of the rescaled variable y = x/ct H for dif- 
ferent Hurst exponents H. The constant c is chosen so 
that y has unit mean. From top to bottom, symbols are 
H = 2/3 (L = 410 3 ), H = 4/7 (L = 1.2 10 4 ) and H = 4/9 
(L = 210 4 ). Dashed lines correspond to the predicted slopes 
of Eq. [3] Dotted line corresponds to H = 1/2, where 
Po(y) = 2y/7rexp(-7ry 2 / 4 )- 



We proceed then to verify that the same scaling expo- 
nent 4> characterizes also the behavior of the conditional 
probability close to the absorbing boundary. This is done 
in Fig. [3l where po(y) is plotted as a function of y for 
some exponents H . For all curves, the rescaled variable 
yo varies in the range 0.1 — 0.4, which a posteriori is ver- 
ified to satisfy the long time limit. Good agreement with 
the proposed conjecture is found. Again, convergence to 
the expected scaling is faster in the superdiffusive case. 
Far from the boundary, po(y) has a Gaussian decay, i.e., 
behaves as the free propagator, as observed above. 

Now, we briefly review some stochastic transport pro- 
cesses where the proposed conjecture is shown to work. 

Levy Flights are Markovian superdiffusive processes 
whose jumps obey a Levy stable (symmetric) law of in- 
dex < n < 2 0, Hoj. Since tt(£) - ^T 1 ^, property 
c) is violated. The Hurst exponent is H = l//i. By 
virtue of the Sparre Andersen theorem [2l|, the persis- 
tence exponent is 6 = 1/2, independent of /i. It has been 
shown that far from the origin po (y) behaves as Levy sta- 
ble law, whereas the conditional distribution close to the 
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hence the exponent 



origin scales as po(y) 
4>= 1/2H = 9/H. 

Within the framework of the so-called Continuous 
Time Random Walk model, a class of subdiffusive walks 
is introduced by assuming that the waiting times between 
consecutive (Gaussian and symmetrical) jumps have a 
power-law decay of the kind <~ 1-Q , with < a < 1, so 
that the Hurst exponent is H = a/2 23]. This process 
is stationary, provided that the walker has made a jump 
exactly at initial time 2J|. In presence of an absorbing 
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FIG. 4: Bounded box of unit size. Conditional probability 
p(x) as a function of x for different Hurst exponents H. From 
top to bottom, symbols are H = 2/3 (L = 410 3 ), H = 4/7 
(L = 1.2 10 4 ) and H = 4/9 (L = 210 4 ). Dashed lines cor- 
respond to the predicted slopes of Eq. [3] Dotted line corre- 
sponds to H = 1/2, where p(x) = tv sin(-7rs;)/2. 

boundary, it has been shown that the survival probabil- 
ity behaves as S(x ,t) ~ x /t a ^ 2 for large times [25|, 
hence the persistence exponent reads 9 = H . Moreover, 
Po(y) vanishes linearly close to the boundary [25|, so that 

4> = i = e/H. 

So far, we have focused on a semi-infinite region; sev- 
eral processes (including translocation) involve however 
finite-size domains [1, [2y| . Asymptotically, for large N 
and close to the origin, the conditional probability in 
a bounded box must be independent of t and xq, i.e., 
Px (x,t) = p{x) and one would expect p(x) ~ x^, where 
<p is the same exponent as the semi-infinite case: the par- 
ticle feels the presence of the other boundary at N only 
when sufficiently close to it. The case of BM (H = 1/2) 
in a bounded box can be analytically solved by address- 
ing the associated Laplacian eigenvalue problem; it turns 
out that for long times 0=1: the conditional probability 
close to the boundaries of the box has the same scaling 
as in a semi-infinite domain. A numerical test for fBm 
with other values of H is provided in Fig. [H which again 
reveals a good agreement between simulation results and 
the proposed conjecture. The same holds true also for 
the other processes considered above. 

Finally, one can also interpret our results via a sim- 
ple heuristic argument. Subdiffusive processes tend to 
explore space more thoroughly than regular diffusion, so 
that close to the boundary their probability of being ab- 
sorbed is larger. We therefore expect <f> > 1 for subdiffu- 
sion. As an example, excluded- volume effects make the 



translocation coordinate subdiffusive, and the untranslo- 
cated polymers are more easily found far from the bound- 
ary at any given time. Analogous arguments lead to the 
conclusion that (f> < 1 for superdiffusion. 
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